Continuous atom laser with Bose-Einstein 
condensates involving three-body interactions 



A. V. Carpentier 1 , H. Michinel 1 , D. N. Olivieri 2 and D. Novoa 1 

1 Area de Optica, Facultade de Ciencias de Ourense, Universidade de Vigo, As 
Lagoas s/n, Ourense, ES-32004 Spain. 

2 Area de Linguaxes e sistemas informaticos, Escola Superior de Enxeherfa 
Informatica, Universidade de Vigo, As Lagoas s/n, Ourense, ES-32004 Spain. 

E-mail: avcarpentier@uvigo.es 

Abstract. We demonstrate, through numerical simulations, the emission of a 
coherent continuous matter wave of constant amplitude from a Bose-Einstein 
Condensate in a shallow optical dipole trap. The process is achieved by spatial control 
of the variations of the scattering length along the trapping axis, including clastic 
three body interactions due to dipole interactions. In our approach, the outcoupling 
mechanism are atomic interactions and thus, the trap remains unaltered. We calculate 
analytically the parameters for the experimental implementation of this CW atom 
laser. 
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1. Introduction 

The experimental realization of Bose-Einstein condensates (BEC) with alkali atoms pQ 
has opened the door to the realization of atomic beam sources of high coherence and 
brightness, with strong analogy to optical laser beams. Such "atom lasers" promise 
unprecedented achievements in fields such as atomic interferometry or gravimetry 
and have been proposed in several configurations, predominantly pulsed [2j [3j H] or 
semicontinuous El [71 El E] , and with different outcoupling mechanisms [TU1 fTTl [T2l [13] . 

For instance, the first atom laser used short radio-frequency pulses to outcouple 
atoms from the cavity, flipping the spins of some of the atoms and releasing them from 
the trap with the aid of gravity [2j. However, despite its success in the first experiments 
on atom laser operation, spin-flipping techniques have some limitations as a versatile 
outcoupling mechanism for coherent matter-wave sources. In first place, there are strong 
fluctuations in the output at high fluxes due to the fact that the spin-flipping mechanism 
populates all accessible Zeeman states |14j . This constitutes a strong drawback for 
practical applications of atom lasers in high precision measurements such as matter 
wave gyroscopes [15]. Moreover, in the former atom lasers, the emission process is 
driven by gravity and thus the resulting matter wave is always emitted in the vertical 
direction. 

Atomic soliton lasers overcome the previous limitations using the mechanism of 
modulational instability to trigger the outcoupling of atoms from the cloud [21 H]- Thus, 
in this case there is no need for gravity since the emission is obtained by an adequate 
combination of nonlinear effects in the atom cloud and the manipulation of the trap. 
In particular, the necessary condition is that the total number of particles in the cloud 
exceeds a critical threshold. Therefore, atomic soliton lasers add arbitrary directionality 
to the emitted beam. However, currently proposed matter-wave sources based on 
nonlinear effects can only operate in pulsed regime, due to modulational instability of 
the beam: a nonlinear effect which yields a burst of atomic solitons [HI [T71 [18] . These 
localized nonlinear waves arise as a result of the perfect balance between dispersive and 
nonlinear effects in one dimensional systems. A set of robust pulses that propagate 
without shape distortion is therefore produced after the emission, thereby precluding 
continuous operation of the matter wave laser. 

The goal of this paper is to present a novel outcoupling mechanism for atomic beam 
sources which overcomes all the previous problems and produces directional continuous 
wave (CW) emission of an atomic beam with constant amplitude. The trick is to use a 
combination of a shallow dipolar potential with attractive two-body and repulsive three- 
body interactions. Such interatomic potentials can be obtained for strongly dipolar 
atomic systems in the presence of adequately tuned magnetic fields. [19] We will show 
that this technique allows both the ability to extract atoms from the trap pE] and to 
overcome modulational instability of the outcoupled beam, thus providing fluent CW 
operation of the matter- wave laser. Moreover, we will show that because of the inclusion 
of the three-body elastic repulsion the emitted wave displays constant amplitude and 
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is very robust against perturbations, permitting the achievement of atomic coherent 
sources with unprecedented brightness and stability. 

2. System configuration and theoretical model 

Let us assume that an elongated BEC is strongly trapped in the transverse directions 
(x, y) by a parabolic potential V± and weakly confined in the longitudinal axis z by a 
shallow dipole trap V z of width L and depth V Q . We also assume a step variation along 
z of the atomic interactions. Thus, for z < (where the condensate is initially placed 
at t = 0) the interatomic forces are set to a negligible value whereas in the z > region, 
the value of the scattering length is such as to yield strong two-body and three-body 
interactions comparable to the trapping force along z. 

We will consider the case of both two-body attractive and three-body repulsive 
forces in the region z > 0. The techniques for controlling the scattering length which 
produces two-body attractive interactions are well-known and may include optical or 
magnetic fields [HJ EH EE]- It has also been proposed that repulsive three-body forces 
can be activated with an optical lattice driven by microwave fields in ultracold gases of 
highly polar molecules [19] . such as RbCs (22] or atoms like Cr (23]. For this technique, 
Hubbard models predict strong nearest-neighbor three-body interactions, whereas the 
two-body terms can be tuned independently. Also, recent experiments with ultracold Cs 
atoms (21] have revealed the existence of the so-called Efimov states, which represent a 
paradigm for universal quantum states in the three-body sector. These remarkable 
results have opened the possibility to novel BEC applications based on three-body 
controlled interactions, as presented in this paper. 

We will consider here a system of iV weakly interacting dipole bosons of mass m, 
trapped in a potential V(r), where the evolution in time t of the condensate wavefunction 
\1/ is correctly described according to experiments (25] by a mean field Gross-Pitaevskii 
equation of the form 

i%— = V 2 ^ + V(r)V + A(z, |*| 2 W, 

at 2m 

where N = J \ ^\ 2 d 3 r. The functions V and A describe respectively the external 
trap and the dipole-dipole induced potential, which can be activated in the region z > 0. 
In this work we consider a cigar-shaped BEC which is tightly confined in (x, y) by means 
of a tubular harmonic trap. A shallow potential provided by an optical dipole trap V z 
[26| [27] keeps the cloud trapped along the z direction. Thus, the explicit form of the 
three-dimensional trap V is given by: 

V(r) = V ± + V z = ^(x 2 + y 2 ) + V (z/L), 

being u± the frequency of the parabolic trap, V the depth of the shallow optical 
dipole potential, and L the characteristic size. The potential barrier can be overcome, 
without destroying it, if there is a spatial variation of atomic interactions. Thus, we 
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will consider the simple case in which the last term of eq. |2] is modulated by a 
step function along z of the form A(z < 0) = and A(z > 0) = f^l^l 2 — f^l^l 4 
indicating U2 = Airh 2 a/m and ^3 = JJib the strength of two-body and three-body 
interactions, respectively. The previous energies are determined by the value of the 
s-wave scattering length a and an effective volume parameter b, which indicates that 
three-body interactions act in a range ~ b 1 ^ 3 . In the present work we will consider a 
and b as positive constants. 

3. One-dimensional eigenstates 

We will consider setups in which the ground state of the optical dipole trap is much 
larger than the ground state of the transverse harmonic potential. In this situation we 
can describe the dynamics of the condensate in the quasi-one dimensional limit as given 
by a factorized wavefunction of the form [T6], [28] \I/(r, t) = &o(x, y) ■ x(z, t). The density 
of atoms per unit length is given by \x\ 2 - Normalizing the squared wavefunction by 
setting \ip\ 2 = y^vralxl 2 , the following dimensionless equation is obtained: 

where the normalized variables in time and space are r = u±t (time measured in units 
of the inverse of the radial trapping frequency) and 77 = z/r± (length expressed in 
units of the transverse size of the cloud r± = ^Jh/mu^) where f(rj) = 2V z /hu±, a 
and f3 = b/Ay2ar\ are the effective trap, two-body and three-body atomic interaction 
coefficients, eq. [3] is known as the cubic-quintic (CQ) nonlinear Schrodinger equation 
due to its dependency with the cubic and quintic powers of the wavefunction. This 
model has been extensively used in nonlinear optics |29j as well as in superfluidity [30] 
and has known analytical solutions for the one- dimensional case. 

In Fig. [1] the geometry of the system to be considered in this paper is shown. The 
shallow potential with depth Vq and width L is represented by the shaded region. Within 
the zone rj > 0, non-zero interatomic interactions can be achieved by magnetically tuning 
the Feshbach resonances [201 CH] or by optical manipulation [211 HH]- 

We have searched numerically for eigenstates of eq. [3] with the form r) = 
<f)(r])e lflT , where [i is the chemical potential. It is well-known that eq. [3]features a cut-off 
in the chemical potential spectrum, /i cr = ^W, which constitutes the upper border of 
the existence domain of localized solitonic solutions [30J. 

Starting from the ground state of the trapping potential [(a) in Fig. [1] and by 
increasing gradually the number of atoms N, the eigenfunctions tend to enter the 
nonlinear zone [(b) in Fig{T] displaying a high spatial localization. Once iV reaches a 
certain critical value the shape of the eigenstate approaches that of a continuous beam 
of fixed amplitude located outside the trap [(c) in FigJT]. In the latter case the value 
of the amplitude A is completely determined by the strength of the atomic interactions 
and can be easily calculated for a plane wave solution of the CQ homogeneous model 
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Figure 1. [Color Online] Amplitude profiles of several eigenstates of eq. [3] calculated 
for different numbers of atoms in the cloud. The right edge of the shallow trap 
(shaded zone) is located at -q = 0. For the zone with 77 > 0, both attractive two- 
body and repulsive three-body interactions are characterized by a = /3 = 1. Inset: 
Dependence of the number of particles in the condensate N on the chemical potential 
fx/fXcr- The dashed (solid) line corresponds to the system displayed in the outer plot 
with a potential of width L [L' = 4L)). Labeled points refer to the eigenfunctions 
displayed. 



As it can be appreciated in the inset of FigHJ in the case of a narrow trap of 
width L, the dashed curve N vs. yu//ii cr monotonically increases (far enough from the 
critical value fi cr where it becomes divergent) indicating that all eigenstates are stable. 
Otherwise, the solid curve corresponding to a system with a trap of width L' = 4L shows 
a region with a negative slope, yielding to the existence of an unstable domain. Notice 
that for /i — > /i cr , both curves merge because the nonlinearity overcomes the trapping 
potential. 



4. Continuous wave emission 



Taking into account the previous eigenstate structure of eq. [3j we have performed 
a systematic analysis of the evolution dynamics of an initial Gaussian atomic cloud 
of unit amplitude and width oj = 5 placed at the center of a shallow square trap of 
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Figure 2. [Color Online] Evolution of a Gaussian atomic cloud in the presence of 
localized gain T = T cr with two-body coefficient a = 1 and different values of the 
three-body atomic interactions: a) = 0. b) /5 = 0.8. c) (3 = 1. In all snapshots, the 
x-axis displays the spatial coordinate rj and the y-axis corresponds to the evolution 
time t G [0, 500]. Below are the final density distribution plots of a) [top], b) [middle] 
and c) [bottom] figures. 

depth Vo = 0.09 and width L = 4, by means of numerical simulations of eq. [3] in 
different parametric regimes. Subsequent feeding of the cloud with more atoms [31] has 
been simulated by adding a linear gain term iTip to eq. [3] (r G 3ft) localized in the 
trapping region. In the absence of nonlinear interactions (a = (3 = 0) the gain effect 
is to increase the density of atoms in the trap. When only two-body interactions are 
considered (/3 = 0) a burst of solitons is emitted towards the region r\ > as predicted 
in previous work pE] [see picture (a) in Fig. [2] . When /3 > and T is set below a certain 
critical threshold T cr which depends on the parameters of the system (in the simulations 
of Fig. [21 we have fixed T = 0.0375 which corresponds to r cr for a = (3 = 1), the first 
soliton emitted becomes wider featuring a flat-top profile [30] which is characteristic of 
systems with CQ nonlinearity [see picture (b) in Fig. Ej. However, if /3 > and T > T cr 
a continuous beam of constant amplitude is outcoupled from the trap [see picture (c) in 
Fig. [2]. Notice that the density plots in Fig. [2] correspond to the final evolution stages 
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(r = 500) of simulations a) [top], b) [middle] and c) [bottom]. 

In terms of the physical description given above, it can be stated that whenever 
T > the system follows the eigenstate structure for growing values of N, but only above 
the critical value T cr the system is able to produce a coherent plane wave of constant 
amplitude A. Furthermore we have checked that the amplitude of the continuous matter 
wave released remains unaltered (unless small scale fluctuations) when three-body losses 
[32] are taken into account. We have numerically verified that whenever losses are 
moderate (typically less than 10 _3 /3) their effect can be compensated by increasing the 
gain in the system. For greater losses the analysis turns out to be more complicated 
and a more careful study, beyond the scope of this paper, must be performed. 




(c) 



(0 



Figure 3. [Color Online] Same as Fig. [5] for different values of the gain coefficient 
r/r cr . In all cases we have fixed a = = 1. In a) T/T cr = 0; b) T/T cr = 0.8; c) 
T/r cr = 0.93; d) T/T cr = 0.98; e) T/T cr = 0.99; f) T/T cr = 1. 



Another important feature of our system consists on the possibility of controlling 
the spatial width of the first soliton outcoupled. For a given pair of nonlinear coefficients 
[a, /3] this can be achieved by varying the gain coefficient T as it can be seen if Fig. [3j In 

a) snapshot of Fig. [3J where T is zero, almost all atoms are emitted in a coherent pulse of 
finite width. The increment in the spatial extension of the atomic soliton becomes non- 
negligible for small variations of the gain parameter in the regime r/r cr rs 1 [see pictures 

b) to f) in Fig. [3]. The case for which the rate of atoms loaded exactly compensates 
losses due to CW emission is illustrated in Fig. 3 — /). It is noticeable that for F > T cr , 
the CW emitted preserves its peak density and constant shape, indicating that the 
method is robust against perturbations. 

In Fig. H] we have plotted the results of many simulations in the whole parameter 
space, showing the working range of the CW atom laser for different a values. Plotted 
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lines define the lower border of the operational region for each a, i.e., line points 
correspond to V — T cr and CW emission is assured for gain values over such limit. 




Figure 4. [Color Online] Working range of the CW atomic laser for three different 
values of the two-body interaction parameter: a — 1.75[black pointed]; a = l[red 
dashed] and a — 0.75[green continuous]. For each value of the three-body interaction 
coefficient (3, it is plotted the corresponding r cr indicating that over this value CW 
emission is always achieved. 



5. Analytic condition for cw emission 

In the following we will derive an analytic expression for the continuous matter wave 
emission based on waveguide modal theory |33j . 

We consider the linear trap (depicted in Fig. [TJas a one-dimensional linear 
step-index waveguide with core index n c = Vq, cladding index n c i = (left trap 
border) and substrate index n s = V n \ (right border), being V n i the effective potential 
created by the nonlinear interactions. For the moment we will treat V n i as a constant 
value. In this context the guiding properties of the structure can be characterized 
by both the normalized frequency f = L^V — V n i and the cut frequency f c (y) = 
arctan^JVni/ (Vo — V n i) + vn, where v is a non-negative integer, which establishes a cut- 
off for the existence of guided modes [33]. The specific case / c (0) corresponds to the 
threshold of the fundamental mode. We assume that if the lowest energy mode is cut 
there are not guided modes supported by the waveguide and just radiation modes can be 
excited. In other words, for f < f c (0) particles cannot be trapped and will eventually 
flow towards r\ > 0. 
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Let us determine the values of V n i for which the non-guiding condition F = 
fo ~~ /c(0) < is fulfilled. Fixing L and Vq to the numbers employed in the simulations, 
we have found the limiting value V£i = 0.05 numerically solving F = 0. Thus for 
V n i > V£i the guide cannot retain the matter wave field, while in the opposite case the 
particles will remain trapped if the proper mode is correctly excited. Hereafter we will 
consider V n i to be field-dependent owing to the existence of nonlinear interactions within 
the region 77 > 0. On the right border of the dipole trap, = a\ijj\ 2 — /Sj^l 4 indicating 
that the potential barrier depends on the particle density. Notice that as \ip\ 2 must 
be real the latter expression predicts a cut-off in the outcoupling mechanism since for 

2 

Vni — fs the emission is forbidden. The agreement between analytical and numerical 
estimations of this cut-off is above 90%. 

We have also determined the dependence of V£i on Vo (fixing L) to be linear with 
slope £ = 0.917, ordinate x — —0.039 and correlation R 2 = 0.9998. Finally, assuming 
that the atomic cloud features a Gaussian transverse profile as in the whole simulations 
we arrive at the analytic expression for the emission condition: 

a + - 4P(ZV + x) 
A > - 2 ; 

where A is the amplitude of the Gaussian atomic cloud loaded in the trap. This 
condition establishes a lower limit for matter wave emission. Continuous operation 
of the atom laser will be assured by keeping the peak particle density A 2 above this 
threshold. In our numerical simulations we have added a linear gain factor in order to 
satisfy this condition. Experimentally this could be realized by continuously loading the 
optical trap with a flux of particles [31] . 

6. Discussion and conclusions 

We have shown in this work that CW operation of a matter-wave laser could be realized 
within an ample parameter space for negative scattering length and repulsive three- 
body elastic interactions. Although the cubic-quintic model can be considered an exotic 
system, it has been proposed that repulsive three-body forces can be activated by 
simultaneously switching on an optical lattice in 77 > containing an ultracold gas 
of particles with high magnetic or electric dipole moment |19j . 

In summary, we have proposed a novel mechanism for continuous operation of a 
coherent matter wave laser. Our system is able to perform a regular and controllable 
emission of coherent atomic beams of constant amplitude with the only limitation of 
the number of particles that can be loaded in the trap. It is important to note that even 
though the present paper is devoted to square traps, our results also apply to other trap 
configurations without loss of generality. We have also shown how the width of the first 
atomic soliton released can be easily controlled by means of variations of the input gain. 
Moreover we have derived an analytic condition for the CW emission which relates the 
peak particle density of the atomic cloud with the tunable parameters of the system. 
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As the techniques for coherently feeding BECs progress our idea could provide a novel 
approach to new types of matter-wave lasers. 
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